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PREFACE 

The  increasing  interest  in  water-based  systems,  such  as  oceanic 
basing  of  missile  and  space-tracking  systems,  stimulates  original 
contributions  of  basic  research  pertaining  to  a  relatively  unknown  and 
relatively  little-understood  operational  environment.  This  Memorandum 
develops  theoretical  models  and  compares  them  with  observed  data  for 


the  Gulf  Stream  and  the  Cromwell  Current. 
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SUMMARY 

Two- layer  and  continuous  models  of  the  Gulf  Stream  and  one- layer 
and  continuous  models  of  the  Cromwell  Current,  respectively,  are  dev¬ 
eloped  from  theoretical  considerations.  Numerical  results  obtained 
from  computer  solutions  of  these  models  are  compared  with  observed  data 
and  are  found  to  show  reasonably  good  agreement.  The  limitations  of 
the  modeling  are  set  forth. 
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LIST  OF  SYMBOLS 

A  *  coefficient  of  horizontal  eddy  velocity 
Ay  =  coefficient  of  vertical  eddy  velocity 

a, a^,a2,a^  =  coefficients 

b, b^,b2,b^  =  coefficients 

b^_  =  known  value  on  the  boundary 
c  =  coefficient 

D,D^,D2  =  thickness  of  a  stream  layer 

Dq  =  thickness  of  the  stream  layer  at  infinity 
D°,D2  =  thickness  of  the  stream  layer  1,  2  at  infinity 

f  =  coriolis  parameter 
g  =  acceleration  of  gravity 

g'  =  g(p2-p1)/p2 
h  =  height  of  the  free  surface 

h^,h2  =  height  of  the  free  surface  of  layer  1,  2 

h  *  dimensionless  variable  h/h 

o 

K  *  constant  -  (pfe  -  pg)/2pb 
K  *  force  vector,  due  to  friction,  etc. 

»v 

L  «  length  of  path 

L  =  as  superscript,  indicates  stream  left  of  axis 

M  =  as  superscript,  indicates  stream  between  stream  L  and 
stream  R 

n  =  ratio  of  the  thickness  of  a  layer  at  the  equator  to  the 

thickness  at  y 
Jo 

o  *  as  subscript  or  superscript,  indicates  at  infinity 
p  =  pressure 

^1*^2  =  Pressure  ^•tl  layer  1»  2 

p  =  as  subscript,  refers  to  "warm  core"  peak 
q  =  7  x  V 

R  ■  radius  of  the  earth 

R  m  as  superscript,  indicates  stream  right  of  axis 
r  =  vertical  co-ordinate  »  (p^  -  p)/(pb  -  ps) 
r(s)  =  vertical  co-ordinate  defining  the  surface 
r(b)  *  vertical  co-ordinate  defining  the  bottom  boundary 


u 

-~r 


-a 

v 


vl’v2 


v  = 


a 

X 


velocity  component  in  x  direction 

unknown  value  in  the  interior 

velocity  with  reference  to  the  earth 

velocity  with  respect  to  the  inertial  system 

velocity  component  in  y  direction 

velocity  components  in  y  direction  of  layers  1,  2 

dimensionless  variable  =  vljl Kgh^ 

horizontal  co-ordinate 

co-ordinate  of  the  "warm  core"  peak 

dimensionless  variable  =  x/y/ 2 


2Kgho/f‘ 


y 


A 

y 


Z(r) 


Zo<r) 


Z' 

o 


horizontal  co-ordinate 
2  1/4 
(2KghQ/B  ) 

critical  fall-off  distance  for  A 

1/4 

dimensionless  variable  =  yy 

function  giving  the  z  co-ordinate  of  water  with  density 
co-ordinate  r 

function  giving  the  z  co-ordinate  infinitely  far  to  the 
right  of  water  with  density  co-ordinate  r 

3Zo(r) 

Sr 


z"  , 

°  ~Ti~ 

dr 

Z'  *  peak  value  of  Z1 
AP 

Z  *  dimensionless  variable  =  Z  /h 
o  o  o 

z  =  vertical  co-ordinate 
z(r)  ■  vertical  co-ordinate,  function  of  r 
z(s)  =  h 


o  =  reciprocal  of  density  (Appendix  B) 
a  =  numerical  coefficient 

3  ■  D^/D^  (two- layer  model  of  Gulf  Stream;  pp.  17-19) 
g  m  (derivative  of  a  coriolis  parameter) 
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Y  “  Ap -^/ AP2  (two- layer  model  of  Gulf  Stream;  p.  18} 

2 

Y  =  0  /g'h0  (one- layer  model  of  Cromwell  Current;  p.  53) 
A  -  z(r)  -  ZQ(r) 

A  =  A/h 

o 

f  dv  du 
^  =  Sx  "  3y 

T]  *  [curl  v  +  7x(u)  x  r)] 

9  =  latitude 

X  =  radius  of  deformation  (Appendix  B) 

_  a2  , 

?  =  y  /z 

p  =  density 

P1,P2,P3  =  density  of  layer  1,  2,  3 
Ps  =  density  at  the  surface 

=  density  at  the  bottom  of  region  of  interest 
p(  )  =  density  function  of 

at  =  density  parameter  =  1000  (p  -  1) 

X  -  q*a7i|( 

y  =  stream  function 

n  *  angular  velocity  of  the  earth 

u)  =  angular  velocity  vector  of  the  earth 
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I.  INTRODUCTION 

k 

This  analysis  applies  the  concept  of  potential  vorticity  to  two 
ocean  currents--the  Gulf  Stream  and  the  Cromwell  Current.  First,  the 
experimental  data  for  the  Gulf  Stream  are  reviewed.  Then  Stommel's^^ 
simple  one-layer  model  of  the  Gulf  Stream  is  discussed,  and  the  assump¬ 
tions  necessary  to  the  theory  are  explicitly  stated.  This  model  con¬ 
tains  no  variable  parameter.  Appendix  A  contains  important  modifica¬ 
tions  that  should  be  scanned  before  a  careful  reading  of  the  body  of 
the  Memorandum. 

GULF  STREAM 

A  two-layer  theory  is  developed  along  the  lines  of  Stommel  which 
gives,  semiquantitatively,  the  features  of  an  inertial  countercurrent 
and  the  recession  of  the  axis  of  maximum  velocity  under  the  Gulf  Stream. 
Note  that  the  theory  is  only  a  consistency  argument  and  does  not  consider 
cause  and  effect. 

Next,  it  is  shown  that  writing  the  equations  of  motion  as  if  the 
stream  headed  due  north  at  the  point  where  the  theory  is  to  be  applied 
gives  almost  no  error. 

The  continuous  theory  is  developed.  The  density  is  used  as  a 
vertical  co-ordinate,  and  in  this  co-ordinate  system  the  assumption  of 
a  level  of  no  motion  leads  to  a  Laplace-type  equation  to  be  solved  in 
a  specified  region  with  given  boundary  values.  It  enables  us  to  cal¬ 
culate  the  density  and  velocity  of  a  cross  section  of  the  Gulf  Stream 

* 

Defined  as  vorticity  of  the  water  divided  by  the  local  gradient 
of  any  conserved  field.  In  this  study,  the  density  is  the  only  such 
field  used. 
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if  we  know  only  the  density  structure  on  the  boundaries  of  the  cross 
section.  First,  a  simple  analytical  model  is  solved  for  a  linear  den¬ 
sity  gradient  above  a  region  of  constant  density.  Next,  a  solution 
that  agrees  well  with  observation  is  obtained  by  computer,  using  ob¬ 
served  data  as  input  data  on  the  boundaries. 

Finally,  a  discussion  of  the  possible  effects  of  friction  on  the 
accuracy  of  the  solution  is  given  in  Appendix  B.  To  discuss  these 
effects,  it  is  necessary  to  derive  the  conservation  of  potential  vor- 

ticity  in  general  form.  The  results  show  that  if  the  Navier-Stokes 

<  5  2  * 

expression  is  used  for  viscous  stresses,  then  for  A^  =  10  cm  /sec 
<  2  ** 

and  A^  =  10  cm  /sec  we  may  get  negligible  effects.  This  completes 
the  investigation  of  the  Gulf  Stream. 


CROMWELL  CURRENT 

.  The  hypothesis  that  the  water  of  the  Cromwell  Current  sinks  at  3° 
north  or  south  latitude,  moves  to  the  equator,  rises  to  the  surface  and 
moves  back  to  the  same  latitude  is  shown  to  be  improbable  because  of 
salinity  data  and  because  of  certain  dynamic  considerations.  It  is 
assumed  that  this  water  comes  from  a  source  extending  across  the  equator 
and  therefore  that  its  potential  vorticity  is  highly  variable.  Simple 
models  are  used  to  show  the  reasonableness  of  this  view.  The  assumption 
is  made  that  the  potential  vorticity  of  a  parcel  of  water  is  proportional 
to  its  distance  from  the  equator.  This  is  because  the  source  of  water 
found  at  some  distance  from  the  equator  is  likely  to  be  farther  from 
the  equator  than  a  parcel  found  on  the  equator. 


Ay  =  Coefficient  of  horizontal  eddy  viscosity. 


** 


Ay  =  Coefficient  of  vertical  eddy  viscosity. 
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Under  this  assumption,  the  continuous  case  is  developed  for  the 
Cromwell  Current,  and  a  computer  solution  gives  the  three  layers  of 
this  case,  one  above  the  other,  semi quantitatively,  using  as  input 
data  the  density  structure  at  the  equator,  on  the  surface  and  at  4° 
latitude.  There  is,  unfortunately,  a  variable  parameter  in  the  theory-- 
which  can  probably  be  measured--to  represent  the  distance  from  the 
equator  of  the  source  of  a  parcel  of  water. 

Adjustment  of  this  parameter  for  best  fit  to  the  data  gives  a 
narrow  source  region  of  width  about  4°  of  latitude.  The  results  are 
rather  insensitive  to  the  value  of  the  parameter.  Finally,  in  Appendix 
B,  the  limitations  of  this  theory  of  the  Cromwell  Current  due  to  fric¬ 
tion  and  vertical  velocities  are  considered.  Here,  in  contrast  to  the 

2 

Gulf  Stream,  values  of  Ay  ■  5  cm  /sec  might  seriously  invalidate  the 
theory. 
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II.  THE  GULF  STREAM- -THE  EXPERIMENTAL  DATA 

A  cross  section  of  the  Gulf  Stream  reveals  that  lines  of  constant 

temperature  and  salinity  rise  as  one  follows  them  from  east  to  west. 

(2) 

This  fact  is  illustrated  in  the  cross  sections  taken  by  Iselin, 
shown  in  Figs.  1  and  2.  In  the  temperature  section  by  Iselin  we  see 
a  typical  feature  of  the  Gulf  Stream,  the  "warm  core"  that  in  this  case 
is  bounded  by  the  20°C  isotherm  near  Station  1227.  Usually,  associated 
with  this  warm  core  is  a  countercurrent  which  is  seen,  for  example,  in 

(3) 

the  velocity  cross  section  computed  by  Neumann  and  shown  in  Fig.  3. 
Neumann's  section  shows  that  the  point  of  maximum  velocity  for  a  given 
depth  recedes  eastward  under  the  Gulf  Stream.  It  is  also  notable  that 
the  maximum  velocity  of  the  countercurrent  recedes  westward.  The  order 
of  magnitude  of  both  these  slopes  is  1/50,  We  also  see  that  both  the 
10°C  and  4°C  isotherm  fall  a  fraction  (1  -  e  *)  of  the  vertical  distance 
between  their  western  and  eastern  elevations  in  a  west-east  distance  of 
about  90  km. 

Figure  4,  by  Iselin,  shows  the  depth  of  the  10°C  isotherm  in  the 
western  North  Atlantic.  These  lines  to  a  large  extent  parallel  the 
actual  transport  lines  of  the  water  feeding  into  the  Gulf  Stream  as 
drawn  by  Iselin.  A  cross  section  of  Iselin's  along  the  42°  meridian 
of  longitude  is  shown  in  Fig.  5. 

THE  ONE-LAYER  MODEL 

Review  of  Previous  Work 

To  illustrate  the  spirit  in  which  this  analysis  will  be  developed, 
let  us  review  briefly  a  one-layer  model  of  the  Gulf  Stream  by  Stomnel. 
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St  o  Hon 


Fig.  1  --  Temperature  section  across  the  Gulf  Stream-- 
Chesapeake  Bay  to  Bermuda 
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Fig.  4  --  Depth  of  the  10°C  isothermal  surface  in  the 

(2) 

western  North  Atlantic  (in  meters)' 


5  —  North-South  temperature  section  through  the  Caribbean 
indicating  approximate  constancy  of  potential  vorticity 
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He  considers  the  cross  section  of  the  stream  as  illustrated  in  Fig.  6. 
There,  h  is  the  height  of  the  free  surface,  D  is  the  thickness  from 
the  surface  to  the  10°C  isotherm,  g  is  the  acceleration  of  gravity, 
p  is  the  density  (gm/cc),  and  f  is  the  coriolis  parameter.  As  a  basis 
for  the  model,  we  assume  the  following  points: 

1.  The  velocity  may  be  considered  to  be  directed  to  the  north, 
and  u  ^  and  v  ^  are  negligible  compared  with  fv  and  g£. 

2.  Frictional  effects  in  the  water  as  it  flows  from  the  Caribbean 
to  the  point  of  observation  are  negligible. 

Assumption  (2)  implies  assumption  (3),  which  is 

3.  The  potential  vorticity  in  a  line  of  constant  density  in  a 
cross  section  is  some  known  function  which  can  be  evaluated 
near  the  source  of  the  stream. 

Assumption  (1)  is  not  serious,  and  the  errors  it  introduces  will 
be  shown  to  be  negligible.  Assumptions  (2)  and  (3)  are  mostly  to  be 
justified  by  the  results,  although  some  discussion  of  them  will  be 
given. 

The  conservation  of  potential  vorticity  is  briefly  derived  in 
Appendix  B.  It  says  that  if  water  flows  from  point  to  point  in  a  layer 
of  constant  density,  and  if  there  are  no  forces,  whose  curl  does  not 
vanish,  acting  on  the  water,  then,  along  streamlines 

—  jr  —  =  constant 


where 


-  dv  du 

c  =  5T  '  S? 

f  =  2Q  sin  0 

Q  =  angular  velocity  of  the  earth 


0  *  latitude 
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Fig. 


6  —  Stommel's  one- layer  Gulf  Stream  model 
(velocity  is  directed  into  the  paper) 
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Stommel's  model  of  the  Gulf  Stream  then  assumes  that  v  and  ^ 

vanish  far  to  the  right  of  the  intersection  of  the  10°C  isotherm  with 

the  surface.  Then,  making  assumption  (3)  in  the  form  of  assuming  that 

the  potential  vorticity  in  the  layer  above  the  10°C  isotherm  is  constant 

in  the  interior,  we  set  it  equal  to  — .  Therefore,  across  the  stream 

o 

we  have 


3v 

iUE  .  f(y) 

D  D^y) 


(1) 


As  to  the  justification  for  assumption  (3)  in  this  particular  case, 
we  first  quote  from  Stommel. 


The  justification  for  applying  the  simple  model 
described  rests  on  the  observed  fact  that  in  the 
central  Atlantic  between  10°  and  35°  N  latitude, 
the  potential  vorticity  of  the  top  isothermal 
layers  (>  10°C)  is  remarkably  constant.  Suppose, 
for  example,  that  the  quantity  DQ  is  the  thickness 
of  the  layer  between  the  sea  surface  and  the  10°C 
isothermal  surface.  The  depth  of  the  10°C  isotherm 
is  shown  in  Fig.  66  [our  Fig.  4],  Since  the  relative 
vorticity  in  central  oceanic  regions  is  small,  the 
reader  may  compute  f/D0  for  the  region  between  10° 
and  35°N  and  convince  himself  that  over  this  large 
part  of  the  ocean  the  potential  vorticity  so  defined 
is  actually  uniform  and  nearly  constant.  '  ' 


Now,  since  this  is  the  region  fron  which  all  the  water  in  the  stream 


is  drawn,  according  to  the  transport  linos  of  Iselin,  it  is  clear  why 
we  may  set  the  potential  vorticity  equal  to  a  constant  across  the  stream 
even  though  it  is  dynamically  a  constant  only  along  stream  lines.''1  Still, 
if  the  water  because  of  turbulent  viscosity  did  not  retain  its  original 
potential  vorticity  as  it  traveled  along  the  stream  lines,  its  value 
might  be  changed  from  its  constant  value  by  the  time  it  arrived  at  the 
observation  point,  and  the  theory  could  not  apply. 
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Next,  we  assume  a  level  of  no  motion  below  the  10°C  isotherm.  So 
we  have,  referring  to  Fig.  6 

(D  -  h)  p2  -  DpL 

The  pressure  gradient  in  the  upper  layer  due  to  the  elevation  of 
the  free  surface  is 


3Pi  ah 

iT  ■  “i  £ 

And,  by  assumption  (1),  all  nonlinear  (u  .  7)u  are  small  and  the 

geos trophic  equation  is  valid. 

,  1  9P1 
fv  =  -  - — 

p  a* 


Combining  the  three  equations  above,  we  get 

,  1  3D 

£v-;  s 


(2) 


where 


g  -  g 


P2  ~  Pi 

P2 


Combining  Eqs.  (1)  and  (2),  we  get 


3  D  1  . 

' 7  (  o> 


(3) 


where 


Imposing  the  condition  that  D  goes  to  Dq  and  v  goes  to  0  at  infinity, 
and  D  ■  0  at  x  =  0,  we  find 

D  =  D  (1  -  e  X/^)  and  v  =  V g  'D^  e  (4) 

o'  *°o 

-3  2 

If  we  set  D  =  800  m  and  g;  »  2  x  10  m/sec  ,  then  we  get  a 
maximum  velocity  of  4  m/sec.  If  we  take  f  at  35°N,  we  get  \  *  48.6  km. 
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Xt  is  of  interest  to  note  here  that  Iselin's  cross  section  gives  X. 
approximately  equal  to  90  km. 

Extending  the  One-Layer  Model--Caribbean  Potential  Vorticity 

It  is  natural  to  try  to  get  more  detailed  models  of  the  Gulf  Stream 
by  elaborating  on  this  simple  model  of  Stommel's.  This  is  what  we  shall 
do  in  the  next  few  pages.  Before  making  these  extensions,  it  is  neces¬ 
sary  to  discuss  assumption  (3)  in  somewhat  more  detail.  We  have  seen 
that  the  potential  vorticity  of  the  water  above  the  10°C  isotherm  in 
the  Caribbean  is  nearly  constant;  this  gave  us  the  justification  for 
the  simple  one-layer  model.  Now,  if  we  are  to  do  many- layered  models, 
or  even  a  continuous  model,  we  must  assume  that  each  layer  considered 
has  constant  potential  vorticity  in  the  Caribbean.  From  Fig.  5  we  can 
see  that  this  is  approximately  true  for  the  water  between  the  24°C  and 
10°C  isotherms.  However,  because  the  thickness  of  the  layers  is  approxi¬ 
mately  constant  below  the  10°C  isotherm,  the  potential  vorticity  increases 
in  this  layer  as  we  move  to  the  north.  But  even  this  will  be  changed 
by  only  a  factor  of  about  two  for  the  source  region  considered,  so  that 
it  still  may  be  profitable  for  qualitative  considerations  to  extend  the 
model  below  the  depth  of  the  10°C  isotherm.  It  is  also  of  interest  to 
note  that  since  we  will  actually  use  density  and  not  temperature  as  our 
marking  property,  the  constancy  of  potential  vorticity  in  the  Caribbean 
is  not  as  good  for  Iselin's  density  cross  sections  as  for  his  temperature 


cross  sections. 
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TWO-IAYER  MODELS 


Establishing  the  Models 

Making  use  of  the  above  comments,  we  shall  now  consider  two-layer 
models  of  the  Gulf  Stream,  as  illustrated  in  Fig.  7.  In  this  model 
(essentially  Stommel's  model  split  in  two)  we  have  two  layers  of  equal 
thickness  of  400  m  at  infinity,  but  of  differing  densities,  and  p 2- 
A  level  of  no  motion  below  D2  gives  the  requirement 

<D2  -  h2>  p3  =  Dlpl  +  D2p2 


and  we  see  that 


hl  + 


where  h^  and  represent  elevations  or  depressions  of  their  respective 
boundaries  from  the  free- equilibrium  height  to  the  west. 

For  the  pressure  variations  we  have 

dpL  Shj^  3p2  dhj^  dh2 

ax-  “  8pi  w  =  8pi  sr  +  g  (p2  pP  air 

Combining  these  slightly  more  complicated  equations  (as  was  done 
in  Stommel's  model)  with  the  approximation 


we  get 


where 


p2  '  P1  p2  P1 


etc. 


'1  f 


2  f 


A 


9D1 

3D2 

Api3T  + 

A  p2  a5T 

8D1 

3D2 

A  p2  + 

A  p2  $T 

pi 

(5a) 

(5b) 


P 


Fig.  7  --  Schematic  of  two-layer  models  of  the  Gulf  Stream 
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We  also  get 


where 


8  A  P2 

— X 


P“Df 


(6a) 


(6b) 


These  equations  hold  wherever  we  have  two  layers.  When  there  is  only 
one  layer,  as  to  the  left  of  x  *  0  in  Fig.  7,  we  need  only  use  the  un¬ 
complicated  one- layer  equations  to  obtain  the  solution  for  that  region. 

It  is  necessary  to  solve  the  system  of  coupled  simultaneous  equa¬ 
tions  and  then  to  impose  boundary  conditions. 

The  solution  of  the  homogeneous  equations  is  of  the  form 


cijX  -cijX  a2*  -a2x 

D.  ■  a.  e  +  a„  e  +  a,  e  +  a,  e 
11  2  3  4 


a,x  -ax  a-x  -a2x 

D-  ■  b,  e  +  b_  e  +  b_  e  +  b,  e 
2  1  2  3  4 

where  the  a's  satisfy  the  deterrainental  requirement 


Y  a 


2 


2 

a 


2 

a 


(7a) 


(7b) 


-0 
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or 


(Y  -  1)  a4  -  \  (Y  +  b  a2  +  -  0 

X2  e  3X 


(8) 


To  determine  the  connections  between  the  eight  unknown  coefficients, 
ai*  bi’  We  su^st^tute  a^ove  solutions  into  either  Eq.  (6a)  or  (6b) 
and  set  the  coefficient  of  each  exponential  equal  to  zero.  For  the 
coefficients  of  the  negative  exponentials,  we  get 


,2  "  al 


2  ‘  a2 


a2  = 


—  b2 


T~b4 


“1  2 

If  we  also  insist  that  the  solution  remain  finite  at  infinity,  we  get 

ax  =  bi  =  a3  =  b3  =  ° 

Then  the  solution  in  the  region  of  two  layers  is,  adding  the  particular 
solutions  D°  and  Dj 


7T  -  a.  -  a0 

.2  1  -a.x  .2  2  -a_x 

_  -  A  1  ,  A.  Z 

D1  =  b2  2  e  +  b4  2  e  +  D1 
a,  a„ 


(9a) 


-a  x  -a  x 

D2  =  b2  e  +  b4  e  +  D° 


(9b) 


If  we  use  the  single- layer  equations  for  the  bottom  layer,  the  solution 
to  the  left  of  x  »  0  is,  if  we  impose  D  =  0  at  x  =  -L. 


(10) 


-x/X  <2li  +  x)A 

.  TN0 

(L^  +  x)/\ 

«  b 

0 

e  -e 

+  D2 

1  -  e 

Now  that  we  have  the  solutions  in  the  various  regions,  we  must  join 


them  by  appropriate  boundary  conditions.  We  shall  always  impose  con¬ 
tinuity  of  layer  thickness  and  of  velocity  at  the  boundaries.  It  should 
be  noted  that  because  potential  vorticity  is  conserved  across  the  stream, 
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^  is  automatically  continuous  across  the  boundaries.  Finally,  by 
comparing  Eqs.  (2)  and  (5b),  we  see  that  the  continuity  of  v  guarantees 
the  smoothness  of  (D^  +  D2)  ,  which  is  the  boundary  line  between  p 
and  p3- 

Examples  of  Two-Layer  Models 

For  our  first  model,  let  us  suppose  that  the  bottom  400-m  layer 

never  rises  to  the  surface  to  the  left  of  x  =  0  and  that  it  goes  to 

its  asymptotic  value  of  400  m.  Then,  if  we  prescribe  as  typical  values 

ApL  =  2  x  10“3  Ap2  =  10*3  D°  =  D°  =  400  m 

-4-1  2 

f  *»  10  sec  g  *=  10  m/sec  X  =  20  km 

we  get  as  a  solution  the  six  parts  of  Eq.  (11) 

D2  “  D2  [l  +  0,236  eXA]  (m) 

D2  «  D°  [l  +  0.618  e-1*618x/^  .  o.382  e~°- 618x/x] 

=  D°  [l  -  0.618  e-°*618x^  .  o.382  e‘1,618x/X]  (11) 

v2  *  0.472  ex/X  (m/sec) 

v2  =  1.236  e-0,618x/x  -0.764  e'1,618x/X 

v^  ■  2.00  e-°,618x/x  +  0.472  e_1*618x/X 
This  is  graphed  in  Fig.  8. 

The  maximum  velocity  is  2.47  m/sec;  there  is  a  velocity  to  the  left 
of  x  *  0  which  was  not  found  in  Stommel's  model,  and  the  maximum  velocity 
in  the  lower  layer  is  displaced  9.62  km  to  the  right.  This  corresponds 
to  a  slope  of  the  axis  of  maximum  velocity  of  between  1/25  to  1/50. 
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Pig.  8  --  Two- layer  model  where  lower  layer  does  not  rise 
to  the  surface,  showing  the  recession  of  the 
axis  of  maximum  velocity  with  depth 
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For  our  next  model,  we  make  the  upper  layer  a  warm  core  because 

we  anticipate  that  this  will  give  a  countercurrent.  However,  it  is 

difficult  to  determine  the  core's  potential  vorticity,  since  none  of 

its  density  lines  extend  to  infinity  on  the  right  where  the  water  is 

at  rest.  Now,  Fig.  9  shows  the  density  structure  which  we  wish  to 

(4) 

model.  The  experimental  section  is  taken  from  Worthington,  his 
Fig.  4. 


From  the  fact  that  X  is  typically  of  the  order  of  40  km,  we  anti¬ 
cipate  that  in  this  model  v  =  0  in  the  middle  of  the  warm  core.  There¬ 
fore,  assuming  that  this  layer  came  from  a  region  of  uniform  potential 
vorticity,  we  may  set  the  potential  vorticity  in  this  region  equal  to 

f/D  ,  where  D  is  the  observed  thickness  of  the  core  at  its  middle 
o  o 

(equal  to  125  m  in  this  case).  If  we  knew  where  this  water  came  from, 

we  could  observe  the  water  in  the  source  region  and  calculate  its 

potential  vorticity.  Then,  as  a  solution,  we  get  Eq.  (12),  which  is 

graphed  in  Fig.  10.  It  is  seen  that  indeed  we  do  get  a  countercurrent 

and  that  the  axis  of  maximum  countercurrent  velocity  recedes  westward 

under  the  stream  as  one  would  anticipate.  The  eight  parts  of  Eq.  (12) 

are  .  . 

x/x  -x/x 

=  39.41  e  -  53.05  e  +440  (m) 


x  -  200 


=  40.48  e 


a 

d”  -  128.6  e 


+  440 


lx  -  200 

1  1  *o  "<*lx/xo 

+  115.5  e 


(12) 


x  -  200 


-88.0  e 


2  X, 


'a2x/xo 

129.0  e  0  +  440 


o-t  =  1000  (/D-l) 


Fig.  9  --  Two- layer  model  to  be  solved  to  give 
countercurrent^^ 


t 


23 


Fig.  10  —  Solution  of  Fig.  9  from  Eq.  (12) 
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-  -102.5  e 


x  -  200] 

“  '  x°  1  -  92.0 


/ 


0-2 

-22.22  e 

x/\. 


x  -  200 


^  /  "a?x/X„ 

-  32.60  e  2  0  +  125 

-x/X, 


v£  -  0.296  e  0  +  0.398  e  0  (m/sec) 


x  -  200 


v*  -  -0.304  e 


■  0.436  e 


al 


x  -  200 


I.  0.388 


(12) 


lx  -  200 


-0.740  e 


a2  X 


'  '  ♦  1.087 


v“  -  -1.190  e 


x  -  200 

X  -a1x/Xc 

1  0  1  +  1.061  e  1  ° 


x  -  200 

®2  X0  -a2x/X0 

-0.880  e  0  ’  +  1.291  e  0 


where 

X  =  28.22  km 
o 

al  -  2.233 
02  ■  0.8934 

This  result  is  only  semiquantitative,  since  Fig.  9  is  a  bad  approxi¬ 
mation  of  the  true  density  structure,  which  has  large  gradients  inside 
the  core. 

ROTATION  OF  CO-ORDINATE  SYSTEM  TO  STREAM  AXES 

Before  moving  on  to  the  continuous  case,  it  would  be  well  to  verify 
whether  assumption  (1)  (page  10)  can  be  valid  in  the  case  of  the  Gulf  Stream. 
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Going  back  to  Stommel's  model,  if  we  include  both  north-south  and  east- 
west  velocities,  we  have,  corresponding  to  Eq.  (1) 

dv  3u 


f  + 


x  ~  f_ 

D  D 


Upon  substituting 


we  get 


f  +  gVf 


fv  -  g 


fu  -  g 


/ 

ax 

,  M 

ay 


-  8  ‘ 


9D 

Sy 


D 


Now,  the  ratio  of  the  third  term  in  the  numerator  of  the  left-hand  side 
of  the  above  equation  to  the  other  two  if  the  stream  is  proceeding  in 
a  north-westerly  direction  is  approximately 


3f 


dD 


5?/f)  57 
1  ills 

f  2 

ay 

which,  by  differentiating  Stommel's  solution  to  get  approximate  values 
for  the  derivatives,  is  of  the  order  where  R  is  the  radius  of  the 
earth. 


Thus,  this  term  may  clearly  be  neglected.  Next,  it  is  well  known 
that  the  Laplacian  equation  is  invariant  under  co-ordinate  rotations. 
Therefore,  in  the  equation  resulting  from  the  above  simplification  we 
may  rotate  the  co-ordinate  system  until  the  derivative  normal  to  the 
stream  contains  all  of  the  numerical  magnitude  of  the  Laplacian  equation. 

If  we  then  call  this  direction  "x,"  we  have  recovered  the  original  equation. 
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THE  CONTINUOUS  CASE 


Development  of  Theory 

This  discussion  follows  in  broad  outline  the  paper  of  Rossby. 
Let  us  define  a  vertical  co-ordinate  r  such  that 


Pb  ‘  P 

r  *  - 

Pb  ‘  Ps 


(13) 


where 

pg  =  density  at  the  surface 

pb  =  density  at  bottom  of  region  of  interest 

Note  that  r  goes  from  0  at  the  bottom  to  1  at  the  surface.  We  also 
define  K,  a  constant,  such  that 


K 


pb  ~  Ps 


2P 


b 


It  follows  that 


p  «  pb(l  -  2  Kr) 


Remembering  Eq.  (1)  for  the  conservation  of  potential  vorticity. 


we  rewrite  it  as 


Sv 

3x 


To  pass  to  the  continuous  case,  we  replace  D  by  an  infinitely  thin  layer 


and  get 


D 


dr 


We  now  define  ZQ(r)  as  a  function,  giving  the  z  co-ordinate, 
infinitely  far  to  the  right,  of  water  with  the  density  co-ordinate  r. 
We  can,  therefore,  rewrite  the  previous  equation  as 


dx 


(14) 


We  should  remember  at  this  point  that  the  above  equation  is  for 


the  conservation  of  potential  vorticity;  as  such. 


dv  .  .  , ,  dv 

g^  is  actually  g^ 


with  the  density  held  constant,  and  so  it  may  be  written 


Equation  (14)  can  be  rigorously  derived  as  an  approximate  form  for 


the  rigorously  derived  Eq.  (65)  in  Appendix  B. 


We  now  define 


A  »  z(r)  -  ZQ(r) 


and  the  equation  above  becomes 


\2H\  _  -  ^ 

dx 


dA 

f-5T 

c 

dr 


Now  the  pressure  at  some  depth  z  is 

_surface 
gPdz 


(15) 


Integrating  by  parts,  we  get 
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where 


p(s)  »  p  on  surface 
p(z)  ■  p  at  z 

Now,  it  is  easily  seen  that 


(H)  ■  (£) ♦  »<■> 

'  'Z  '  '0  \  / O 


(17) 


To  obtain  the  first  term  on  the  right-hand  side  of  the  above  equa¬ 
tion,  we  substitute  the  derivative  of  the  right-hand  side  of  Eq.  (16). 
Equation  (17)  becomes 

p(s) 


(I)  -  «<•>  (§)  •  «<■>(*£)  -(b)  f ’< 

v  '2  '  'P  x  'p  x  'P  J 

P(z) 


gzdp  +  gp(z) 


(18) 


The  second  and  fourth  terms  cancel,  and  the  third  term  becomes 

p(s) 


-gz(s) 


8p(s) 


/ 

p(z) 


r  Sz(p) 

8  -§5T"  dp 


Now,  since 


this  becomes 


:(s)  -  h,  m  ||  and  dp  ,  .2Kpbdr 


2Kpbg  i 


■N*)  ♦  fi 


r(s) 

© 

BA 


dr 


So,  finally 


r(s) 

BA 

55T 


dr 


(19) 


Note  that  Eq.  (19)  is  also  correct  if  we  replace  x  with  y.  This 
will  be  of  use  later  in  this  Memorandum  in  an  analysis  of  the  Cromwell 


Current. 
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At  this  point,  it  is  convenient  to  introduce  a  level  of  no  motion 
at  some  arbitrary  great  depth.  For  convenience  we  take  it  to  be  the 
greatest  depth  attained  by  If  this  level  is  taken  to  be  the  zero 

point  for  the  z-axis,  then  we  have 
z(r(s)) 


■ 


gpdz  +  Pfcgz  (r(b))  =  constant 


z(r(b)) 


(20) 


Integrating  by  parts  we  get 


0  -  p(s)  gz(r(s))  -  p(b)  gz(r(b))  -  J  gzdp  +  p(b)gz(r(b)) 

r(b) 


r  (s) 

-  J°&zdp 


Canceling  the  second  and  fourth  terms  and  then  taking  I  5-  )  gives 


U1 


p(s)g 


Sz(r(s)) 

3x 


r(s) 

2K*ho|z(r(.))5|M  +  f  gdrj 

r(b) 


dr  >  =  0 


Substituting  this  in  Eq.  (19)  cancels  the  unwanted  first  and  second 
terms,  and  adding  the  integrals  gives 


(!l),  •  /i dr 

r(b) 

and,  finally,  the  geostrophic  relationship  gives 

r(z) 


(21) 


(22) 


This,  then,  is  the  formula  which  we  may  use  to  calculate  the  vel¬ 
ocity  when  we  have  calculated  the  final  A(x,t). 
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Taking  the  derivative  of  Eq.  (22)  we  find 


- 

.  2K£  j 

UxL 

*  J 

s2a 

— y  dr 
dx 


(22a) 


r(b) 


and  then 


a_ 

dr 


2  Kg  3  A 

'  '  i7 


(23) 


But  if  we  differentiate,  with  respect  to  r,  Eq.  (14)  for  the 
conservation  of  vorticity,  we  may  set  the  right-hand  side  equal  to  the 
right-hand  side  of  Eq.  (23)  ,  and  upon  rearrangement  of  terms  there 
results  finally,  with  primes  indicating  differentiation  with  respect  to  r 

3A 


2«*Z'o  ih 


2" 

o 


dr 


dx 


M 

d* 


(24) 


This  is  the  final  equation.  It  now  becomes  necessary  to  define 
the  region  in  which  it  is  to  be  solved.  Since  r  is  our  vertical  co¬ 
ordinate,  the  region  is  bounded  below  by  r  =  0.  On  the  surface  far  to 
the  right  r  =  1,  but  as  we  move  left,  we  see  from  Fig.  1  that  the  lines 
of  density  come  to  the  surface  and  thus  the  surface  value  of  r  must 
decrease.  Therefore,  the  region  in  which  Eq.  (24)  is  to  be  solved  must 
look  something  like  Fig.  11. 

It  should  be  remembered  here  that  while  the  v's  in  Eqs.  (23)  and 
(14)  are  both  in  the  geometrical  co-ordinate  system,  Eq.  (14)  is  only 
an  approximation  to  Eq.  (65)  in  Appendix  B,  and  therefore  Eq.  (24)  is 
also  an  approximation. 
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Fig.  11  --  Boundary  of  the  region  to  be  solved  in  (x,r) 
co-ordinates  for  Eq.  (24) 
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Boundary  Conditions 

Next,  we  consider  the  boundary  conditions  which  must  be  satisfied. 
Clearly,  to  the  far  right  the  boundary  condition  must  be  A  **  z(r)  -  zQ(r) 
=  0.  Now,  far  to  the  left,  if  potential  vorticity  has  been  conserved 
and  the  velocity  vanishes,  the  profile  of  density  can  only  have  been 
displaced  upward,  and  so  A  =  constant.  In  between,  on  the  free  surface, 
we  see  that  where  the  density  r  reaches  the  surface,  the  boundary  con¬ 
dition  states  that  A  must  equal  the  distance  which  that  density  line 
has  risen.  That  is,  A  =  ZQ(1)  *  zQ(r)  =  w*iere  h  is  the  height  of 
the  free  surface.  But  the  height  of  the  free  surface  is  negligible 
compared  to  zc(l)  -  zQ(r)>  and  therefore  the  surface  boundary  condition 

is  A  =  ZQ(1)  “  z0(^)* 

There  remains  only  the  bottom  boundary  condition.  From  Eq.  (22) 
we  see  that  the  velocity  for  r  =  r(b)  is  zero.  This  is  true  all  along 
the  line  of  constant  density,  r  =  r(b).  From  Eq.  (22a)  we  also  see 


p*p(b) 


0 


But  this  implies  by  Eq.  (14)  that  =  0.  Therefore,  the  boundary 
condition  on  the  bottom  is 


3A  =  o 


(25) 


For  purposes  of  computation  it  is  convenient  to  introduce  dimension¬ 
less  variables 

A  X 


/\ 
V  ' 


1y2Kghc 


2  -  Z  /h 
o  o  o 


(26) 
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In  terms  of  these  new  variables,  Eqs.  (22)  and  (24)  become 

r 

A  ~~  , 

v  =  -  |  t—  dr 


..A 

eA 

-J 

^  1 

r(b> 

A 

* 

L" 

O 

’  T> 

o 

(27) 

A 

dr2 

The  nondimensional  problem  to  be  solved  is  sketched  in  Fig.  11. 

It  is  to  be  noted  that  with  the  exception  of  the  core,  the  data  necessary 
for  solution  comprise  only  the  density  gradients  on  the  surface  to  give 
the  surface  profile  and  boundary  condition,  and  the  density  in  vertical 
sections  at  infinity  to  the  right  to  give  Z'  and  Z 


Solution  for  Continuous  Case 

Analytical  Model.  It  is  of  interest  now  to  obtain  the  simplest 
possible  analytical  solution  for  the  continuous  case.  We  consider  a 
linear  density  gradient  in  which  the  entire  field  of  variable  density 
comes  to  the  surface  at  one  point.  To  see  how  we  shall  represent  this, 
consider  the  problem  given  by  the  dotted  line  in  Fig.  12.  Here  the  lines 
of  constant  density  come  to  the  surface  linearly  over  a  small,  finite 
distance.  Along  this  line,  A  is  as  seen  in  the  projection  on  the  left. 

If  now  we  pass  to  the  limit  where  the  small  finite  distance  goes  to  a 
point,  we  may  consider  that  the  A  projection  remains  unchanged.  In  the 
interior  of  the  resulting  region,  signified  by  the  solid  line,  we  have 
a  simple  Laplace  equation,  and  on  the  boundary,  simple  conditions.  This 
problem  is  easily  solved,  and  the  solution  is 


A 

A 


00 

r'  i 

r (2n-l)nr 

/  2  cos 

-i  (2n-l)Z 

CM 

_ 1 

-(2n-l)TTx 

2 


e' 


(28) 
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Fig.  12  --  Solution  in  (x,r)  co-ordinates  for  a 
linear-density-gradient  Gulf  Stream 
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Through  Eq.  (27)  the  velocity  is 


A  plot  of  this  solution  is  seen  in  Fig.  13,  where  the  lines  of 
constant  density  and  velocity  are  drawn  and  also  a  line  showing  the  axis 
of  maximum  velocity  with  depth.  The  mean  slope  of  the  axis  of  maximum 
velocity  is,  in  dimensionless  co-ordinates,  about  1/1.27.  The  theory 
naturally  fails  at  the  surface  at  x  =  0  because  of  the  infinite  density 
gradient  and  at  the  bottom  because  we  have  placed  the  level  of  no  motion 
at  too  shallow  a  depth  to  keep  from  straining  the  linear-density-gradient 
approximation.  The  failure  in  the  graph  of  the  bottom  dotted  line  to 
come  to  the  surface  shows  the  effect  of  plotting  only  three  terms  of 
the  series.  The  ratio  of  the  vertical  scale  to  the  horizontal  scale  is 

D 

o 

We  shall  use  somewhat  more  precise  constants  to  extract  the  most 

possible  from  the  theory.  From  Fig.  1  we  pick 

D  =  900  m 
o 

f  =  0.81  x  10"4  sec-1 

=  0.8  x  10"3 
P 

We  find  the  slope  of  maximum  velocity  to  be  1/46.5. 

Since  the  maximum  nondimens ional  velocity  is  0.743,  we  find  the 
maximum  velocity  to  be  1.99  m/sec.  Finally,  we  find  K  ■  33.1  km.  A 
plot  of  the  nondimensional  surface  velocity  is  also  seen  in  Fig.  13. 
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Finally,  to  extract  the  greatest  possible  amount  of  information 
from  Eq.  (24),  a  numerical  solution  was  made.  The  left-hand  edge  was 
taken  as  Station  1231  in  Fig.  1  and  the  right-hand  edge  taken  near  Sta¬ 
tion  1223.  The  region  in  which  Eq.  (24)  is  to  be  solved  in  (x,r)  co¬ 
ordinates  is  shown  in  Fig.  14.  Note  the  large  "warm-core"  peak  on  the 
top  of  which  we  find  r  values  of  1.83.  This  is  to  be  compared  with  the 
value  of  1.00  for  r  at  infinity  on  the  right.  In  this  core,  just  as  in 
the  core  of  the  two-layer  model  of  the  Gulf  Stream  previously  solved,  it 
is  necessary  to  specify  the  potential  vorticity.  This  is  done  rather 
arbitrarily  by  assuming  that  in  the  core,  x—  *  0  directly  under  the  peak. 
Thus,  this  gives  the  potential  vorticity  in  the  core  as  f/Z',  where  Z^  = 

Z'  is  just  at  the  peak.  Notice  that  this  assumption  is  contradicted  by 

Sv 

the  results  of  the  computation  which  give  ^  #  0.  We  could  undoubtedly 
arrive  at  a  consistent  set  of  values  in  the  core  by  iteration,  but  they 
would  not  alter  the  qualitative  results  except  by  smoothing  out  the  grad¬ 
ients  near  the  peak.  The  vertical  line  under  the  peak  is  now  regarded  as 
another  line  of  origin  just  like  a  vertical  line  at  infinity  in  the  east, 
from  which  we  measure  A  for  the  water  in  the  core.  Now  it  is  clear  that 
adding  a  constant  number  to  all  A  above  some  given  r  cannot  change  the 
velocity  structure.  It  is  equally  clear  that  in  the  numerical  solution 
it  will  be  necessary  to  have  A  be  a  continuous  function  in  the  interior. 
Therefore,  all  the  mesh  points  beneath  the  highest  point  in  the  core 
must  and  may  have  the  same  value  of  A  as  the  line  of  greatest  r  which 
manages  to  pass  beneath  the  peak  and  extend  to  infinity  on  the  right. 

If  this  value  of  r,  forming  the  bottom  of  the  core,  is  rQ,  then  the 

constant  value  of  A  mentioned  above  is  seen  to  be  Z(x  ,r  )  -  Z  (r  ) 

p  o  o  o 

■  A  ,  where  x  is  the  x  co-ordinate  of  the  peak.  This  value  of  A  must 
op  o 
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<x 


Final  Gulf-Stream  problem  to  be  solved  in  (x,r)  co-ordinates 


-39- 


then  be  given  to  each  point  beneath  the  peak.  To  obtain  the  surface 
boundary  condition  for  densities  within  the  core,  we  add  to  this  constant 
value  the  distance  which  that  density  line  must  come  up  to  reach  the 
surface.  Note  that  the  many  arbitrary  elements  in  this  treatment  of 
the  core  are  noted  in  the  two- layer  core  solution,  page  21. 

The  density  structure  in  Iselin's  section  at  Station  1223  on  the 
east  is  shown  in  Fig.  15.  This  density  structure  was  approximated  by 
the  three  straight  lines  as  shown.  These  straight  lines  gave  the  value 
of  Z^  except  at  the  intersections.  If  a  mesh  point  in  the  numerical 
solution  is  close  to  such  an  intersection,  then  the  value  of  Z'  is  taken 
as  the  mean  of  the  two  on  either  side.  Also,  in  the  equation  at  such 
a  point,  the  third  term  of  Eq.  (24)  comes  into  play.  To  evaluate  Z", 
the  difference  of  the  Z'q  above  and  below  were  divided  by  what  seemed, 
on  inspection  of  the  actual  graph  of  density  versus  depth,  to  be  the 
proper  value  of  6r.  To  obtain  Z'q  in  the  core,  the  graph  of  density  in 
Iselin's  Station  1228  down  to  200  m  as  shown  in  Fig.  16  was  used.  By 
the  above  procedure  we  obtained  the  coefficients  of  Eq.  (24),  the  region 
in  which  it  is  to  be  solved,  and  the  boundary  conditions. 

Suppose  that  we  construct  a  mesh  for  the  numerical  solution  of  the 
problem  and  that  there  are  a  total  of  n  mesh  points  including  the  inter¬ 
ior  and  the  boundary.  If  we  then  express  Eq.  (24)  in  finite  difference 


form,  then  we  have  n  equations  of  the  form 


S  a.  u<  =  Z  c,  br 
j  XJ  J  r  £r 

l  =  1>  2,  ...,  n 


where  Uj  represents  the  unknown  values  in  the  interior, 


values  on  the  boundary;  the  a's  and  c's  are  known  coefficients.  Since 


i 
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Fig.  15  --  Density  structure  at  Station  1223  In  Fig.  1 
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Fig.  16  --  Density  structure  at  Station  1228  In  Fig.  1 
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the  a's,  b's,  and  c's  are  all  known,  then  for  each  l  the  right-hand  side 
is  just  a  number,  and  we  may  write  Eq.  (29)  as 


£  %  “i 


1 )  2  j  •••)!! 


(30) 


u  =  d 


So  if  the  matrix  could  be  inverted,  we  would  have  the  solution. 

Now,  if  the  numbering  of  the  mesh  points  is  done  in  a  regular  manner, 

=> 

it  is  easy  to  see  that  in  the  matrix  A,  all  the  terms  will  be  zero  except 

1/2 

those  in  a  band  whose  width  is  of  the  order  N  around  the  main  diagonal. 
Computer  Solution.  If  we  devise  a  computer  program  which  takes 

1/2 

advantage  of  these  zeros,  a  400  x  400  matrix  may  be  stored  in  400  x  400 

2 

=  8000  cells  instead  of  (400)  =  160,000.  In  such  a  case  one  need  perform 

no  calculations  with  the  zeros  off  the  diagonal  band,  enabling  fast 
computation.  Such  a  program,  using  Gaussian  elimination  and  written 
by  J.  F.  Holt  of  Space  Technology  Laboratories,  solved  the  mesh  indicated 
in  Fig.  13  on  an  IBM  7090  in  one  minute. 

In  Fig.  14  we  have  indicated  the  values  of  Z'  and  Z"  derived  from 
the  density  data  of  Figs.  15  and  16.  The  surface  profile  comes  from 
Iselin's  surface  densities.  The  large  peak  represents  the  warm  core. 
First,  a  solution  was  attempted  in  which  the  solid  line  represented  the 
boundary  of  the  region.  This  solid  line  is  the  exact  transcription  of 
Iselin's  surface  data.  The  result  of  those  calculations--that  there  is 
a  sharp  decrease  in  A  at  the  extreme  left  edge- -ref lects  the  fact  shown 
in  Fig.  1  that  isotherms  18°C  through  12°C  rise  near  the  surface  around 
Station  1229  but  do  not  reach  the  surface  until  Station  1231.  It  is  of 
interest  that  this  side  lobe  of  core  water  is  not  present  in  all  of 
Iselin's  section,  where  frequently  the  density  lines  come  directly  to 
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the  surface.  We  might  then  conjecture  that  the  lobe  is  an  anomaly  not 
reflecting  the  true  dynamic  structure  of  the  stream.  The  dotted  line, 
therefore,  is  drawn  to  indicate  roughly  what  the  boundaries  would  look 
like  if  the  isotherms  mentioned  above  came  directly  to  the  surface. 

In  support  of  this  elimination  of  the  core  lobe,  note  that  in  the 
solution  performed  with  the  solid  boundary,  almost  all  the  points  between 
the  dashed  boundary  and  the  solid  boundary  had,  for  a  given  r,  larger 
values  of  A  than  the  boundary  value  for  that  r.  That  is  to  say,  the 
solution  forced  these  lines  of  density  to  rise  above  the  free  surface-- 
which  is,  of  course,  physically  impossible.  The  author  feels  that  this 
is  the  correct  solution  to  the  equation  and  does  not  represent  a  computa¬ 
tional  error  or  a  result  of  using  finite  differences. 

In  the  solution  (retaining  the  core),  as  mentioned,  the  lines  of 
constant  density  (which  largely  parallel  the  lines  of  constant  tempera¬ 
ture)  slope  steeply  to  the  right  of  Station  1231  due  to  the  small  boundary 
values  of  A  near  Station  1231.  The  true  situation,  seen  in  Fig.  1,  is 
that  the  density  lines  are  rather  flat  in  this  region.  This  certainly 
suggests  that  we  should  throw  out  the  lobe  and,  shifting  the  solid  line 
to  the  right  as  a  first  guess,  perform  the  solution  with  the  dashed 
lines  as  boundary.  Values  for  parameters  used  in  the  solution  were 

h  =  2000  m 
o 

2K  =  1.79  x  10"3 
f  *  0.81  x  10  4  sec  ^ 

It  is  the  solution  performed  within  the  dashed- line  boundary  which 
is  written  in  schematically  at  a  few  of  the  total  of  400  points  in 
Fig.  14.  This  solution  has  been  transformed  to  (x,z)  co-ordinates  in 
Fig.  17  which  may  be  favorably  compared  with  Fig.  1,  the  experimental  data. 


200 


Fig. 


17  —  Solution 
dotted  line 


in  (x,z)  co-ordinates 
in  Fig.  14  (Z^(r)  per 


of  the  region  bounded 

«*•  5  *  0) 


by  the 
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However,  there  are  still  a  few  delicate  points  in  this  solution. 

First,  the  inconsistency  of  the  assumptions  with  the  results  in  the  core 
reflects  our  ignorance  about  the  core's  source  region.  The  inconsistency 
could  undoubtedly  be  corrected,  simply  if  tediously,  by  the  iterative 
procedure  mentioned  before.  This  would  not  change  the  qualitative 
solutions;  it  would  in  fact  be  only  a  formal  manipulation.  Secondly, 
we  find  that  two  points  next  to  the  bend  in  the  dotted  line  rise  just 
slightly  out  of  the  water.  However,  as  shown  in  Fig.  17,  the  error  is 
small  and  in  context  can  only  mean  that  we  did  not  quite  guess  the  true 
dynamic  situation.  With  reference  to  this  point,  the  whole  pattern  in 
Fig.  17  shows  a  shift  somewhat  to  the  left  with  reference  to  Fig.  1. 

This  could  in  large  measure  be  removed  if  the  obtuse  angle  in  the  dashed 
line  in  Fig.  14  were  moved  right  to  form  a  right  angle,  showing  the  lines 
of  constant  density  coming  to  the  surface  in  a  very  small  region.  This 
would  also  alleviate  the  problem  of  their  coming  out  of  the  water. 

Since  the  density  lines  in  Fig.  17  follow  the  observations  so  closely, 
there  is  little  question  that  the  velocities  will  agree  with  the  observed 
ones,  since  we  determine  the  velocities  in  this  theory  from  the  density 
field  by  directly  applying  the  geostrophic  equation. 

Interestingly,  the  distance  that  the  lines  of  constant  density  fall 
to  their  equilibrium  height  is  accurately  given,  as  is  the  slope  of  the 
axis  of  maximum  current  under  the  stream.  We  see  that  the  core  introduces 
a  "plateau"  in  the  deeper  isotherms  even  though  in  the  (x,z)  co-ordinates 
the  core  looks  rather  insignificant. 

One  final  discrepancy  in  the  results  should  be  mentioned.  The 
Isotherms  do  not  rise  to  the  proper  height  on  the  left.  This  is  simply 
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due  to  a  slightly  erroneous  boundary  condition  inserted  into  the  computer. 
It  could  be  easily  remedied,  but  the  agreement  is  sufficient  as  it  stands 
for  one  to  be  assured  of  the  theory. 
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III .  CROMWELL  CURRENT 


Figure  18  shows  a  density  cross  section  of  the  Cromwell  Current. 
Figure  19  shows  a  salinity  cross  section  at  the  same  longitude  in  May 
of  1958. 

THEORY  OF  THE  CROMWELL  CURRENT 

We  will  try  to  construct  a  theory  of  the  Cromwell  Current  using 
ideas  of  potential  vorticity.  An  earlier  attempt  in  this  direction  was 
made  by  Fofonoff  and  Montgomery. They  suggested  that  water  flowed 
toward  the  equator  in  a  layer  of  constant  thickness,  conserving  potential 
vorticity;  then  at  the  equator  it  rose  to  the  surface  and  spread  to  the 
north  again  only  to  sink  and  again  move  back  toward  the  equator. 

Now,  Fig.  19  shows  that  there  are  widely  varying  salinities  in  a 
cross  section.  These  data  do  not  suggest  that  there  is  a  cross-stream 
flow;  they  suggest  jets  unconnected  laterally  and  moving  east  and  west. 

/  g\ 

This  is  verified  in  the  report  of  King,  which  shows  similar  salinity 
cross  sections  at  140°,  120°,  and  110°W.  Further,  if  the  northern  and 
southern  waters  mixed  at  the  equator,  then  the  concept  of  potential 
vorticity  of  a  water  mass  in  this  sinking  and  rising  water  would  be  even 
more  nebulous  than  it  is  now;  we  might  also  expect  to  find  uniform  sal¬ 
inity  throughout  the  section,  since  it  would  all  have  passed  through 
the  mixing  process.  On  the  other  hand,  if  the  water  did  not  mix,  and 

*  (6) 

The  author  is  indebted  to  John  Knauss'  for  this  drawing  and  for 

all  other  density  data  on  the  Cromwell  Current  mentioned  in  his  article. 

The  data  were  obtained  on  the  H.  M.  Smith  at  140°W  longitude,  April  17 

to  24,  1958. 
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Latitude 

Fig.  18  —  Density  cross  section  of  the  Cromvjell  Current 


Fig.  19  --  Salinity  cross  section  of  the  Cromwell  Current 
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yet  this  circulation  took  place,  then  we  might  expect  to  see  a  variation 

in  time  of  the  salinity  values  at  the  equator  as  different  parcels  of 

water  moved  by.  Such  variations  are  absent  in  the  data  of  repeated 

(6)* 

stations  at  the  equator  by  Knauss. 

It  is  true  that  the  water  at  the  equator  is  an  approximate  mean  in 
salinity  between  the  northern  and  southern  water.  Although  this  could 
show  mixing,  it  indicates  that  the  water  in  the  stream  comes  from  a  slow- 
moving  area  in  the  west,  where  there  is  a  similar  regular  variation  of 
salinity  across  the  equator. 

It  is  of  interest  that  this  salinity  cross  section,  as  drawn  by 

Knauss, does  not  have  the  same  extremely  regular  structure  as  the 

oxygen.  Graphs  showing  phosphorus,  in  addition  to  oxygen  and  salinity, 

(8) 

are  also  among  those  of  King,  which  are  reproduced  by  Knauss.  This 
may  show  that  the  oxygen  and  phosphorus  distributions  are  governed  not 
by  vertical  mixing  but  by  the  Cromwell  Current's  effect  on  some  biological 
processes  which  do  not  affect  salinity. 

To  gain  a  rough  dynamical  insight  into  whether  or  not  it  is  reason¬ 
able,  in  view  of  the  newer  data,  to  say  that  water  moves  in  the  manner 
described  by  Fofonoff  and  Montgomery,  we  shall  construct  a  simple  model. 

We  note  first  from  Fig.  18  that  in  the  region  of  swiftest  current  the 
isopycnics  are  practically  linear  between  0°  and  3°N  latitude.  So  we 
suppose  that  a  layer's  thickness,  represented  by  h,  is  given  by 

h  =  (n  -  1)  hQ  (1  -  y/yQ)  +  hQ  (31) 

* 

Despite  these  arguments,  recent  measurements  by  Knauss  (personal 
communication)  indicate  such  large  (~  20  m/sec)  north-south  velocities 
in  the  Cromwell  Current,  supporting  the  arguments  of  Fofonoff  and 
Montgomery. 


-51- 


where  h  is  the  thickness  at  y  =  y  ,  and  nh  is  the  thickness  at  y  =  0. 

We  see  that  n  is  the  ratio  of  the  thickness  of  a  layer  at  the  equator 

to  its  thickness  at  y  .  We  shall  select  y  =  3°N  in  this  case.  The 

o  o 

equation  for  the  conservation  of  potential  vorticity  then  is 


(32) 


We  note  immediately  that  if  we  apply  this  equation  to  the  south 
side  also,  replacing  yQ  by  -y  ,  then  at  the  equator  where  y  =  0  there 


is  a  large  unrealistic  discontinuity  in 


i*u 

ay* 


Inserting  Eq.  (31)  into  Eq.  (32)  and  integrating  with  the  boundary 
condition  u  =  0  at  y  =  y  ,  we  get 


is 

2 


u  =  -s-  (y  -  y  ) 


(33) 


The  maximum  value  of  n  in  the  central  core  as  picked  off  Fig.  18 
is  about  3.  We  then  find  that  the  maximum  expected  velocity  is  about 
400  cm/sec.  This  is  too  large.  Therefore,  arguments  against  the  north 
south  flow  of  water  are  the  discontinuity  of  ~  and,  secondarily,  the 
large  velocity  resulting  from  this  assumption. 

We  might  then  suppose,  in  harmony  with  the  salinity  considerations 
that  the  water  came  from  a  region  where  the  water  was  moving  slowly  and 
the  isotherms  were  horizontal  across  the  equator.  In  this  case,  the 
water  would  come  from  a  region  of  varying  potential  vorticity.  For 
example,  the  water  at  the  equator  would  have  zero  potential  vorticity, 
in  agreement  with  observation. 

To  describe  this  we  would  write  in  place  of  Eq.  (32) 

3u 


9y~aV  By 
h  h 

o 


(34) 
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Carrying  through  the  same  manipulation  as  in  Eqs.  (31),  (32),  and 
(33),  we  find 

<35> 

This  time  the  maximum  velocity  is  about  90  cm/sec,  which  is  closer 
to  observation,  and  there  is  no  discontinuity  at  the  equator.  Further¬ 
more,  when  n  <  1,  as  actually  occurs  below  200  m,  we  get  a  countercurrent. 
The  model  breaks  down  near  the  surface,  since  the  lines  of  constant  den¬ 
sity  are  no  longer  straight  lines.  However,  we  see  that  n  <  1  (roughly) 
there  also. 

Now,  it  should  be  noted  that  if  we  write  the  equation  of  potential 
vorticity  in  the  form 


a£i 

h 


(36) 


o 

we  cannot,  for  reasons  of  continuity,  expect  the  equation  to  be  correct 

with  a  =  1  because  the  water  must  come  in  from  the  side  as  the  stream 

moves  from  west  to  east,  with  its  layers  thickening  and  velocities 

increasing;  a  then  represents  the  latitudinal  extent  of  the  source  region. 

o 

For  example,  if  a  =  2,  we  will  imagine  that  the  water  found  at  2  N 
began  its  rapid  movement  toward  the  equator  at  4°N.  We  remember  in 
this  regard  that  the  source  region  of  the  Gulf  Stream  was  measured  in 
thousands  of  kilometers,  while  the  width  of  the  stream  itself  was  measured 
in  hundreds  of  kilometers. 

If  we  carry  through  the  model  above  with  an  a,  we  find 
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Here,  if  we  put  <j  *  1.25,  we  find  the  maximum  velocity  to  be 
140  cm/sec.  This  small  value  of  a  suggests  that  the  stream  brings  in 
water  from  near  the  equator  over  a  long  longitudinal  distance. 


ONE -LAYER  MODEL 

We  now  construct  a  one- layer  Cromwell  Current  model  analogous  to 
Stommel ' s  Gulf  Stream  model.  We  use  Eq.  (36)  with  the  geostrophic  equa¬ 
tion 


fJyu 


& 
8  9y 


and  get 


ay2  y  ay  8  ho 

Introducing  as  dimensionless  variables 


J.  2 

v- 


(37) 


(38) 


h  » 
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yy 


1/4 


b2 

V  8  * 

g  ho 


we  get 


A  fr2h  a  ah  a4  a  a4 
v  ■*  ~  -  v  ^  -  uy  h  =  -  y 


y72  -  y  9 
ay 


The  solution  to  this  is 


(39) 


i  -  .  e-W2  ♦  b  X"  +  ,  (40> 

Discarding  the  growing  exponential  for  finiteness  at  infinity,  we  get 

1/2  a2 


b  e 


•a 


7/2 


+  1 


A  A 

Setting  h  »  n  at  y  »  0,  we  find  b  ■  n  -  1;  therefore 


(41) 


1/2  a2/7 

h  -  h  (n  -  1)  e"a  y  /2  +  h 
o  o 

, -  1/2  a2 

u  -  a  ^  (n  -  1)  *JT\  e'a  :! 


(42) 
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These  results  suggest  substituting  §  *  y  /2  in  Eq.  (39).  This  gives 

A 

-  ah  =  -1  (43) 

3? 

We  see  that  we  have  recovered  the  sort  of  equation  we  solved  for 
the  Gulf  Stream;  thus  we  will  plan  to  make  this  substitution  in  the 
continuous  case  also. 


Continuous  Case  of  the  Cromwell  Current 

The  equation  corresponding  to  Eq.  (15)  is 

3u  Pya  If 

=  (a.  -  Dpy 


(44) 


As  we  noted  previously,  Eq.  (19)  holds  also  in  the  y  direction, 
so  we  have,  also  analogous  to  Eq.  (21) 


( 

\dy/; 


z  "pfu  ”  "2Kgpb  /  fy  dr 


A  A 


(45) 


So  in  this  case 


.  2Kg 

0y 


K 
J  3y 


dr 


(46) 


Despite  the  additional  term  in  Eq.  (44),  the  analysis  follows  as  before. 
Differentiating  the  right-hand  side  of  Eq.  (45),  first  with  respect  to 
y  and  then  with  respect  to  r,  and  equating  it  to  the  right-hand  side  of 
Eq.  (44),  differentiated  with  respect  to  r,  we  get 


(9) 

*Robinson  has  justified  the  use  of  the  geostrophic  relation  in 
the  second  momentum  equation  within  2  to  10  km  of  the  equator. 
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Making  the  substitutions 


(_jL) 

\  2Kgh  / 


-1/4 


o 

A  z 

7  0 

Z  =  — 
o  h 

o 


we  get,  using  primes  for  r  differentiation 

JL  2lA  .LM,  „  I1 

2  ~  a2  a3  a 

y  dy  y 

a2 

Under  the  substitution  §  *  y  /2  this  becomes 

A 


2*  z>  Z" 

+  _2  .  °  M  _ 

7  +  -  ’  F  ar  - 

o 


A 

7  »  A 

O  £A 

A,2  dr 


St  S5 

The  equation  for  the  velocity  becomes 


r 

“  "  V21C8h  o  J 


£« 


(48) 


(49) 


And  we  see  that  we  have  precisely  the  same  mathematical  problem  as  for 
the  Gulf  Stream. 


Cromwell  Current  Data 

In  the  case  of  the  Cromwell  Current  it  will  be  necessary  to  specify 
the  surface  boundary  condition,  A  =  0  at  infinity  to  the  north,  and  A(r) 
at  the  equator.  To  specify  A(r),  two  accurate  density  sections  are 
needed — one  at  the  equator  and  another  around  4°N  or  S.  For  the  section 
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at  the  equator,  five  stations  from  the  R.  V.  Horizon,  April  7  to  11, 

1958,  were  averaged,  and  the  result  is  plotted  in  Fig.  20.  For  the 

station  to  the  north  from  which  A's  are  to  be  measured,  four  H.  M.  Smith 

stations  from  3°  and  3.5°N  and  S  latitudes  were  plotted  on  one  graph, 

and  a  smooth  line  was  drawn  between  them.  This  is  seen  in  Fig.  21. 

To  obtain  A  as  a  function  of  p>  a  value  of  p  was  selected  in  Fig.  21 

and  also  found  in  Fig.  20.  The  difference  in  depth  was  divided  by  1000  m, 
A  A 

thus  giving  A(r).  In  Fig.  22,  A  is  plotted  as  a  function  of  r  at  the 
equator. 

In  Fig.  23,  we  see  a  plot  of  the  region  in  which  we  are  to  solve 

Eq.  (49)  in  (r,0)  co-ordinates,  together  with  the  surface  profile  drawn 

from  pertinent  H.  M.  Smith  stations.  On  the  left  side,  the  solid  line 
A  A ,  A 

shows  A(r);  Z^  and  Zq  are  also  seen  as  computed  from  the  straight  line 
drawn  in  Fig.  21.  In  the  mesh  which  is  to  be  solved,  there  are  285 
boundary  and  interior  points. 

Before  proceeding  to  the  complete  solution,  it  might  be  profitable 
to  make  a  simple  continuous  model.  For  this  purpose  we  model  A(r)  at 
the  equator  by  the  cosine  function  shown  by  the  dotted  line  at  the  left 
in  Fig.  23.  In  this  simple  model  we  imagine  a  linear  density  gradient 
and  also  a  surface  boundary  condition  to  be  imposed  by  the  solution. 

It  is  easily  seen  that  an  appropriate  solution  satisfying  the  above 
conditions  is 


A  ■  A  cos  2nr  e 
o 


-2n?a 


1/2 


(50) 


u  -  Kgah  Aq  sin(2nr)  e  2rr^a 


1/2 
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Fig.  20  —  Average  of  five  stations  0°N  Lat. ,  140°W  Long. 


28 


Depth  ( m) 


CTt  =  1000  (/>H) 


59 


Fig.  22  —  A  at  0  N,  140  W--based  on  Figs.  20  and  21 
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Fig.  23  --  Region  for  solution  of  Cromwell  Current  in  (r,0)  co-ordinates 
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We  see  that  the  characteristic  fall-off  distance  for  A,  y  is 
1/2 

obtained  by  setting  2rr5a  -  1,  or 


where 


ych  "  1/2 

na 


r  2Kgh 


1/4 


To  keep  from  straining  the  linear -density-gradient  approximation 
more  than  we  have  already,  suppose  we  take 


h  =  700  m 
o 

2Rg  =  3.6  x  10* 
8  »  2.3  x  10 


•13 


-1  -1 
m  sec 


Then 


y  ■  456  km 
Jo 


258 

*ch  *  1/4 

a 


km 


Now  the  observed  value  for  y  ^  is  about  150  km.  This  suggests  that 
a  might  be  about  9.  But  this  neglects  an  important  fact  which  we  can 
see  in  Fig.  23.  For  most  values  of  r,  Z^  is  about  1/20  instead  of  1, 
which  was  implicit  in  the  assumption  of  a  linear-density  gradient.  The 
smaller  value  of  Z^  will  make  the  constant-density  lines  fall  off  much 
more  rapidly  than  before.  We  see  that  a  linear-density  approximation 
is  extremely  poor  and  cannot  be  trusted.  Similar  conclusions  are  obtained 
by  evaluating  the  velocities  due  to  this  simple  model.  However,  one 
point  worth  noting  is  that  there  is  no  westward  surface  current  in  this 
model.  This  is  due  to  a  mistreatment  of  the  surface  boundary  condition 
in  this  simple  model,  which  assumes  that  the  surface  has  constant  density. 
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Computer  Solution 

Now  we  continue  with  the  complete  solution  of  this  model  of  the 
Cromwell  Current.  The  problem  is  as  sketched  in  Fig.  23.  The  actual 
co-ordinates  in  which  Eq.  (49)  is  solved  numerically  are  (§,r)  co-ordi¬ 
nates.  This  requires  spacing  which  was  much  more  variable  in  the  %  co¬ 
ordinate  than  that  used  in  the  x  co-ordinate  for  the  Gulf  Stream  if  we 
hope  to  get  values  for  the  current  both  near  and  far  away  from  the  equator, 
using  a  reasonable  number  of  mesh  points.  For  example,  0.1°  corresponds 

—  Lf  o  “1 

to  §  =  2  x  10  ,  and  5  corresponds  to  5  *  5  x  10  .  As  a  result  of 

the  highly  variable  spacing,  the  elements  of  the  matrix  for  inversion 

varied  by  a  factor  of  10^.  However,  extra  precision  was  used  in  the 

inversion,  the  results  are  smooth,  and  there  is  no  reason  to  suspect 

errors.  Perhaps  it  should  be  noted  here  that  a  priori  there  is  no  good 

reason  to  expect  excellent  results  from  this  theory,  even  with  a  as  an 

adjustable  parameter  because,  for  example,  there  is  no  reason  to  suppose 

that  all  the  water  at  a  given  latitude  in  the  stream  originated  at  the 

same  latitude.  Two  facts  save  us:  a,  in  agreement  with  the  discussion 

after  Eq.  (36),  is  close  to  1;  and  the  solution  is  relatively  insensitive 

to  the  value  of  a-  In  fact,  as  can  be  seen  from  the  simple  model  just 

1/4 

presented,  the  decay  distance  of  A  is  roughly  proportional  to  a  ,  and 

1/2 

the  velocity  to  a 

The  region  in  which  the  equation  was  solved  extended  to  §  =  2,  or 
about  10.15°.  In  the  numerical  solution,  all  the  values  of  A  had  declined 
by  factors  of  at  least  20  by  5°,  so  that  the  solution  effectively  did 
extend  to  infinity  away  from  the  equator. 

In  Fig.  24  the  complete  solution  is  drawn  in  (Z,9)  co-ordinates. 

By  comparing  with  Figs.  25  and  26  from  Knauss,^  we  see  that  the 
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Latitude 


Pig.  24  -  Numerical  solution  of  Fig.  23  expressed  in  (Z,0)  co-ordinates 


Depth  (m) 
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Latitude 


7  6  5  6  4  3  2 


Station 

Fig.  25  —  Observed  velocity 
cross  section  at  140°W^ 
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Velocity  (cm/sec) 


Fig.  26  —  A  composite  east-west  current 
profile  at  the  equator  at  140°W^ 
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undercurrent  is  displaced  somewhat  upwards  from  observation,  that  the 
velocity  is  too  slow,  and  that  the  decay  distance  is  slightly  too  large. 
These  last  two  effects  can  undoubtedly  be  corrected  by  introducing  an 
a  of  about  1.4.  In  agreement  with  our  first  simple  dynamical  model, 
this  shows  that  the  source  of  the  water  in  the  Cromwell  Current  is  not 
very  far  from  the  equator,  and  thus  that  the  current  must  grow  slowly 
as  it  moves  across  the  Pacific--not  all  at  once  in  a  few  degrees  of 
longitude. 

The  upward  displacement  of  the  undercurrent  may  be  an  effect  of 


vertical  friction,  as  discussed  in  Appendix  B.  Or,  since  the  density 

lines  so  closely  parallel  the  true  ones,  it  might  show  a  failure  of  the 

geostrophic  equation,  =  ^  However,  it  is  barely  conceivable  that 

it  could  also  be  due  to  a  combination  of  small  errors  in  the  calculation 
A 

of  A  at  the  equator  (Fig.  22),  poor  values  of  r  on  the  surface  between 
0°  and  0.5°,  and  finally  a  true  level  of  transition  between  the  under¬ 


current  and  the  deep  current  at  300  m  instead  of  400  m  as  drawn  by 


Knauss.  That  this  is  possible  we  see  from  the  error  limits  on  the  deep 


velocity  as  drawn  in  Fig.  26. 

Perhaps  it  should  be  emphasized  at  this  point  that  there  is  no 

A 

boundary  condition  in  this  problem  forcing  A  to  go  to  0  anywhere  short 
o  A 

of  10  N,  and  so  the  fact  that  A  goes  to  0  more  or  less  as  observed  is 


a  nonimposed  result  of  the  theory, 

Hie  theory  gives  larger  values  than  observed  for  a  surface  current 
between  1.5°  and  2.5°.  This  is  a  direct  reflection  of  the  large  surface 


density  gradient  and  is  extremely  sensitive  to  the  surface  density  data. 
Nevertheless,  there  is  not  much  room  for  manipulation  of  the  data 


-67- 


available  at  the  present  to  escape  this  result.  Therefore,  these  large 
values  are  a  real  result  of  the  theory.  However,  since  the  calculated 
density  lines  follow  so  closely  the  observed  ones,  it  seems  difficult 
to  explain  the  discrepancy  except  by  a  failure  of  the  geostrophic  equation. 
Perhaps  they  reflect  an  Ekman  layer.  These  large  values  of  the  velocity 
would  only  be  increased  by  increasing  a* 

One  explanation  might  be  that,  as  we  see  in  Fig.  18,  some  of  the 
density  lines  rise  around  3°  as  they  move  away  from  the  equator.  This 
behavior  has  been  excluded  by  the  theory  and  yet  might  slow  down  the 
surface  current.  The  effect  may  also  show  the  influence  of  the  wind 
near  the  surface. 

The  sporadic  behavior  of  the  surface  current  near  the  equator 
possibly  results  from  a  lack  of  surface-density  data,  to  which  the 
surface  current  is  extremely  sensitive.  A  linear  extrapolation  between 
0.5°  and  0°  was  used,  and  it  gave  these  strange  results.  Perhaps  better 
surface-density  data  would  give  the  proper  currents. 

In  agreement  with  observation,  some  density  lines  around  100  m 
reach  a  maximum  elevation  around  1°  and  then  decline  toward  the  equator. 

The  density  lines  below  400  m  are  about  5  m  from  their  asymptotic 
elevation  at  x  =  ®  when  they  run  off  Fig.  24. 
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IV.  CONCLUDING  REMARKS 


Excellent  agreement  with  observation  has  been  obtained  using 
potential- vorticity  theorems  with  no  arbitrary  parameters  for  both  the 
Gulf  Stream  and  the  Cromwell  Current.  These  results  indicate  small 
values  for  the  eddy-viscosity  coefficients.  It  is  predicted  that  the 
Cromwell  Current  grows  slowly  as  it  crosses  the  Pacific  and  does  not 
appear  full-blown  in  a  few  degrees  of  longitude,  nor  does  the  water  in 
it  move  toward  and  away  from  the  equator  in  the  manner  suggested  by 
Fofonoff  and  Montgomery.  This  should  have  important  implications  for 
a  causal  theory.  The  problems  of  the  detailed  dynamics  of  these  streams 
have  therefore  been  shown  to  be  equivalent  to  predicting  the  density 
structures  on  the  surface  and  in  one  or  two  sections,  and  to  knowing 
some  characteristics  of  the  source  region.  There  is  some  need  for  more 
detailed  data  at  the  surface  of  the  Cromwell  Current  to  improve  the 
calculations . 

Even  though  the  theory  is  noncausal,  it  could  provide  a  framework 
for  other  studies.  For  example,  our  simple  continuous  model  of  the 
Gulf  Stream  might  be  useful  for  an  analytical  study  of  internal  waves 
on  the  Gulf  Stream,  or  in  a  stability  analysis. 
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Appendix  A 

THE  CONTINUOUS-MODEL  CORE 

Stommel's  model  (pp.  4-14)  can  be  extended  to  a  fully  three- 
dimensional  model  by  at  least  two  equivalent  analytical  techniques. 

One  is  to  satisfy  the  Bernoulli  integral  on  the  stream- line  which  runs 
along  the  coast.  Another  is  to  insist  that  the  flux  moving  into  the 
coast  below  a  given  latitude  equal  the  flux  moving  north  across  it  in 
the  Gulf  Stream.  Either  requirement  determines  the  depth  of  the  upper 
layer  along  the  coast  and  thus  a  complete  solution.  When  the  layer 
thickness  goes  to  zero  we  have  the  cross-section  of  Eq.  (4).  In  models 
of  more  than  one  layer,  satisfying  the  Bernoulli  integral  will  give  the 
relationship  at  a  given  latitude  between  the  boundary  thicknesses  of 
the  several  layers.  In  this  Memorandum  we  have  assumed  these  values 
to  be  known  from  observation.  It  would  be  a  coincidence  if  we  picked 
values  which  gave,  for  any  latitude,  a  proper  cross  section  for  a  solu¬ 
tion  which  satisfied  the  Bernoulli  integral.  In  fact,  it  can  be  shown 
analytically  that  the  models  on  pages  19  and  33  do  not  satisfy  the 
Bernoulli  integral  at  any  latitude.  However,  they  still  illustrate  the 
types  and  quality  of  results  obtainable  by  the  use  of  constant  potential 
vorticity  and  are  valuable  on  that  account.  Also,  the  continuous  cal¬ 
culation  in  Fig.  17  and  Fig.  27  so  precisely  uses  observed  boundary 
conditions  that  one  may  hope  that  detailed  calculation  would  show  the 
Bernoulli  integral  to  be  satisfied.  Thus,  in  the  main,  this  work  must 
be  regarded  as  a  preliminary  exploration  of  the  types  of  results  perhaps 
obtainable  in  more  rigorous  fashion  by  dynamical  boundary  conditions. 

The  good  qualitative  agreement  obtained  gives  encouragement  to  further 
efforts  to  construct  such  a  theory. 


3 


29  28  27  26  25  24  23 

Station 

27  —  Solution  in  (x,r)  co-ordinates  of  the  region  bounded  by  the 

dotted  line  in  Fig.  14  (5  times  Z'(r)  per  Fig.  14;  L  0) 

o  ox 
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The  treatment  of  the  core  in  the  continuous  model  on  pages  37-39 
is  unnecessarily  forced  and  arbitrary.  It  has  been  replaced  by  another 
treatment,  and  the  continuous  calculation  repeated.  The  result  is  shown 
in  Fig.  27. 

To  discuss  this  new  technique,  let  us  set  the  vertical  scale  by 
assuming  Zq(1)  =  0.  (See  p.  27) 

Then  r 


Z 

o 


(r) 


/ 


Z'(r)  dr 


From  the  definition  of  A  before  Eq.  (15),  we  see  that  at  the  surface 

the  boundary  condition  is  A  =  -ZQ(r).  (Remember  that  z^r)  =  0  at  the 

upper  boundary.)  Thus  we  may  determine  the  boundary  condition  on  the 

water  in  the  core  from  the  observed  values  of  Z^(r)  in  the  core. 

When  the  values  of  Z^(r)  shown  in  Fig.  16  were  used  in  calculation, 

the  resulting  solution  was  physically  impossible.  The  solution  was 

physically  reasonable  only  for  value  of  Z^(r)  in  the  core  from  4  to  6 

times  that  shown  in  Fig.  16.  A  factor  of  5  was  used  for  the  calculations 

leading  to  Fig.  27.  This  is  consistant  with  the  facts  that  the  water 

in  the  core  is  advected  from  lower  latitudes  and  that  there  is  a  large 

horizontal  velocity  shear  in  the  core,  of  such  a  sign  that  the  observed 
dz  SZo 

is  smaller  than  .  Since  we  cannot  precisely  identify  the  source 
of  the  core  water  there  is  no  reliable  way  of  evaluating  this  multipli¬ 
cative  constant  within  the  present  theory  except  to  say  that  it  must  be 

greater  than  unity  for  -r-  <  0.  (See  Eq.  (1)).  The  fact  that  the  present 

ox 

calculation  gives  such  a  factor  is  encouraging. 
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Appendix  B 

DERIVATION  OF  POTENTIAL- VORTICITY  EQUATION 

This  appendix  gives  a  derivation  of  the  potential-vorticity  equation 
after  the  manner  of  Elliason  and  Kleinschmidt^*^  and  also  gives  esti¬ 
mates  of  the  error  due  to  friction. 

DERIVATION  OF  POTENTIAL-VORTICITY  EQUATION 

The  equation  of  motion  in  an  absolute  inertial  system  is 

dv 

JZ  =  -  7  (p  -  a  7  P  +  K  (51) 

dt  ~  ^  ~ 

where  p  is  the  pressure,  a  the  reciprocal  of  the  density,  7cp  the  forces 

derivable  from  a  potential,  and  K  all  the  other  forces,  such  as  friction. 

Equation  (51)  may  be  written 

dv  2 

+  (7  x  v)  x  v  +  7(1/2  v  )  *  -  7cp  -  o.7£  +  K  (52) 

Taking  the  curl  of  Eq.  (52)  we  get 

aq 

■—  +  7  x  (q  x  v)  =  (7  x  K)  +  S  (53) 

Ot  ^  ~  ~  ~ 

where  q  =  7  x  v  and  S  =  7a. x  7p.  Equation  (53)  may  be  expanded,  and 

~  ~ 

by  using  7  ■  (7  x  v)  and  recombining  we  get 

4?  +  q  (7  •  v)  -  (q  x  7)  v  »  (7  x  K)  +  S  (54) 

If  we  define 

X  =  q  *  a  7  t  (55) 

then 
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We  have  the  identity 


d<zi>  a* 

—  -  7(jl)  -  W  •  TV 


(57) 


Note  that  the  last  term  on  the  right,  like  the  others,  is  a  vector. 
Expanding  ^  by  Eq.  (54)  and  substituting  Eq.  (57),  Eq.  (56)  becomes 

,d+ 


Since  -;a  *  a(7  •  v)  is  the  equation  of  continuity, 
d  t 


(58) 


^  *  aVt  •  < (S  +  7xK)  +  (q  •  7)v 
at  I  ~  ~  ~ 


q  -|a(7(^)  -  7V  .  7*)|  (59) 


It  may  be  easily  verified  that 

7 ♦  •  (q  •  7)v  =  q  •  (7v  •  7*) 

If  ♦  is  conserved,  then  ^  =  0;  and  if  ♦  =  g(p),  then 


Since 


then 


And  if 


7t  •  S  S  0 


S  =  7(~)  x  7p 

p  — 


^  -  a  7*  •  (7x  K) 


(60) 


(61) 


then 


♦  “  P 


4r  (curl  v  •  —  7p)  =  —  7p  •  (7  x  K) 
dt  ~  p  o  ~ 


(62) 


And  if 


then 


dE  (curl  Z  *  p  *£> 


(63) 
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Now  in  fact  we  assert  ♦  =  p  and  K  =  0.  So  we  may  define  a  co¬ 
ordinate  system  whose  "vertical"  co-ordinate  is  this  density.  Define 
the  rate  of  change  of  this  co-ordinate  for  a  given  particle  of  water 
as  its  "vertical"  velocity.  This  is  clearly  zero  if  the  density  is 
conserved.  Now  define,  in  the  layers  of  constant  density,  vectors 
perpendicular  to  the  gradient  of  the  density.  This  constructs  an  orthog¬ 
onal  co-ordinate  system.  In  this  co-ordinate  system,  7p  has  only  one 
nonzero  component,  which  is  in  the  same  direction  as  the  vertical  com¬ 
ponent  of  curl  v.  Since  we  have  an  orthogonal  system,  the  dot  product 
in  Eq.  (62)  gives  only  one  term,  the  product  of  the  "verticai"  component 
of  the  curl  and  the  gradient  of  the  density.  Perhaps  we  might  emphasize 
at  this  point  that  the  vertical  discussed  above  is  not  generally  in  the 
same  direction  as  the  gradient  of  the  geopotential,  or  the  vertical 
defined  by  the  radial  component  of  the  coriolis  force.  Nevertheless, 
this  is  a  completely  general  theorem.  The  only  approximations  are  that 
friction  may  be  neglected,  that  density  surfaces  may  be  defined,  and 
that  density  is  conserved  along  streamlines. 

Returning  to  Eq.  (62),  curl  v  is  the  absolute  vorticity  with  respect 

to  an  inertial  system.  Let  us  call  the  velocity  with  reference  to  the 

Earth  u,  and  that  with  respect  to  the  inertial  system  ,  for  the  pur- 

pose  of  this  brief  discussion.  Then 

u  =  (u  +  is  x  r)  (64) 

~a  — r  ~ 

where  ou  is  the  Earth's  angular  velocity  and  r  is  the  Earth's  radius. 

So 

(7x  u  )  =  (7  x  u  )  +  7  x  (<jo  x  r)  (65) 

Inserting  this  into  Eq.  (62)  we  see  that  7p  picks  out  the  component  of 
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the  coriolis  force  parallel  to  itself.  The  value  of  this  component  is 
of  course  very  close  to  sin  0. 

If  we  call  this  component  f  and  perform  the  dot  product,  then 
Eq.  (63)  becomes 


d_ 

dt 


|"(curlH  v  +  f) 


(^Vp) 


} 


0 


(66) 


where  curl,,  v  is  the  curl  of  the  horizontal  velocity. 

In  this  Memorandum  we  apply  the  theorem  to  specific  cases  as 

follows:  In  the  source  region  we  observe  the  value  of  —  .Jcu-irl  u 

0  i  1  ~r 

+  7x(crx 


Now  in  any  of  the  co-ordinate  systems  mentioned  above,  this  dot 
product  has  one  term--the  product  of  the  "vertical"  components--which 
is  much  larger  than  the  others.  So  no  matter  which  of  these  systems 
we  pick,  we  get  essentially  the  same  number  for  the  observed  potential 
vorticity  if  we  take  only  the  product  of  the  vertical  components. 

Now  let  us  imagine  that  we  follow  the  water  flowing  and  we  apply 
the  theorem  in  the  correct  co-ordinate,  system  to  guarantee  that  when 
we  look  at  7p  •  (curl  u  +  7  x  u>  x  r)  again  at  some  new  point  along  the 
streamline  we  will  find  the  same  value.  Then  at  the  new  point  we  shall 
again  write  approximate  equations. 

All  the  actual  mathematics  in  the  body  of  this  Memorandum  deals  with 
these  approximate  equations;  there  is  no  need  to  explicitly  treat  the 
highly  distorted,  time-varying-density  co-ordinate  system.  Nonetheless, 
it  should  be  kept  in  mind  that  it  is  the  conservation  theorem  holding 
in  this  co-ordinate  system  which  enables  us  to  write  the  approximate 
equations.  (There  may  be  cases  in  geophysics,  and  especially  in  model 
experiments,  where  the  difference  in  verticals  is  not  negligible.) 
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Note  especially  that  using  the  theorem  in  this  manner  makes  it  independent 
of  the  hydrostatic  assumption. 

For  an  example  of  application  we  treat  the  Gulf  Stream  system.  In 
its  source  region  we  observe  that  the  distance  between  two  isotherms, 
divided  by  the  local  coriolis  parameter,  is  about  a  constant.  In  evalu¬ 
ating  this  constant,  it  makes  no  practical  difference  which  of  the  three 
possible  verticals  we  use;  so  we  write  simply  f/(dz/dp)  =  constant,  along 
lines  of  constant  density.  Now  we  imagine  that  the  correct  theorem 
"carries"  the  water  to  some  cross  section  of  the  Gulf  Stream.  And  so 
we  write,  along  a  line  of  constant  density, 


dp 


(67) 


In  the  first  pages  of  this  Memorandum,  potential  vorticity  is 
conserved  in  layers  up  to  800  m  thick.  Note  that  if  we  integrate  Eq.  (67) 
between  two  densities  and  consider  the  velocities  to  be  constant  in  the 
range  of  integration,  we  get 


where  D  is  the  thickness  of  the  layer. 


Errors 


The  equation  for  conservation  of  potential  vorticity,  including 


friction,  is 


.  f  i  I  d(T|  •  -  7p)  . 

dt  \  p  7p  *  ^CUrl  v  +  7x(uj  x  r)]  >  =  — ~  dg  ~  ~  7p  •  (7  x  £)  (69) 


where  T1  =  [curl  v  +  7x(<o  x  r)] 

Using  Eq.  (69)  we  may  now  try  to  estimate  the  possible  effects  of  friction 
on  the  foregoing  theory.  We  estimate  when  the  cumulative  effect  of 
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(7  x  K) ,  acting  over  the  path  of  the  water  between  the  source  and  the 
point  of  observation,  would  be  large  enough  to  invalidate  the  theory. 
As  a  first  approximation  we  might  multiply  T|  by  the  time  necessary  for 
the  water  to  travel  this  path. 

GULF  STREAM 

Let  us  consider  vertical  friction  in  the  Gulf  Stream.  We  have 


Therefore, 
is  L/v 


an  -  n^A 


a3y 

v  3xdz2 


(70) 


we  have,  if  the  time  necessary  for  travel  from  the  Caribbean 


21 

7 


A  x  10 
v 


(71) 


where 


L  -  length  of  path  ■  2000  km 
Dq  =  depth  =  500  m 
v  ■  velocity  =  10  cm/sec 
X  *  radius  of  deformation  =  50  km 

From  this  we  see  that,  for  the  values  of  Av  currently  seen  in  the  liter- 

2 

ature  (e.g.,  5  or  10  cm  /sec),  the  friction  effect  is  small. 

For  the  horizontal  eddy  viscosity 

3 

T)  ~  K,  3  (72) 


This  implies 


i  X2v 


Ah  x  10 


(73) 
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So  we  see  that  large  values  of  the  horizontal  eddy  viscosity  would 

5  2 

invalidate  the  theory,  and  we  must  have  A  <  10  cm  /sec  at  least. 


CROMWELL  CURRENT  ERRORS 

Typical  parameters  for  the  Cromwell  Current  are 


D  =  200  m 
L  =  104  km 
u  =  30  cm/sec 


(74) 


A.  =  100  km 

Then  for  the  vertical  eddy  viscosity 

1  D2u 
o 

2 

So  we  see  that  a  vertical  coefficient  of  5  or  10  cm  /sec  might  be 
serious,  especially  in  the  more  slowly  moving  portions. 

For  the  horizontal  eddy  viscosities  we  obtain  much  the  same  conclu¬ 
sion  as  for  the  Gulf  Stream  (i.e.,  A,  <  10^),  since 


^4^3^10^  (76) 

'  X  u 

Of  course  it  should  be  emphasized  that  these  order-of-magnitude  estimates 
do  not  conclusively  show  that  any  effect  is  unimportant  or  important 
because,  for  example,  the  expression  for  friction  is  quite  hypothetical, 
and  any  given  effect  might  be  highly  damped,  or  indeed  be  nonlinearly 
coupled  and  therefore  produce  a  large  effect. 
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